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1. INTRODUCTION 
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For more than thirty years, quantum chromodynamics [QCD] has been established as the theory of the strong 
interaction, yet fundamental questions still remain unanswered. Foremost amongst these are the issues of confinement 
' and the details of the hadron spectrum. As was recognized early on, studies in Coulomb gauge (where the framework 
is naturally connected to physical degrees of freedom) should provide for a promising arena of endeavor. However, 
\^ ■ what was also recognized at the time was that the inherent noncovariance of Coulomb gauge made technical progress 
£^ \ difficult This is one reason why the linear covariant gauges (such as Landau and Feynman gauges) have been 
. preferred for both perturbative and nonperturbative calculations. Clearly though, since the physical world is gauge 
invariant, it is certainly worthwhile to study the different gauges to gain more insight. Perturbation theory is one area 

1 ^ ', where results in different gauges can be unambiguously compared. In addition, perturbative results form the basis of 

many nonperturbative studies since they provide a reliable way of dealing with ultraviolet divergences. 

In the last few years, Coulomb gauge studies have been enjoying significant progress. There exists an appealing 
picture for confinement: the Gribov-Zwanziger scenario 0, y, 0] whereby the temporal component of the gluon 
\ propagator provides for a long-range confining force whilst the transverse spatial components are suppressed in the 
infrared (and therefore do not appear as asymptotic states). Recent lattice studies seem to confirm this Q. A 
Hamiltonian-based approach [1, 0, H, 0, E3] also exists and proves adept at describing various physical features of 
the system [ill, EH, EH • More pertinent to this study, results in the Lagrange-based (Dyson-Schwinger) functional 
' integral a ppr oach (which is especially suitable for perturbative calculations) to Yang-Mills theory have recently become 
OO ! available [ll EE El, [13, El] • 

In this work, we consider the quark contributions in the first order functional approach to Coulomb gauge QCD, 
using results obtained in the pure Yang-Mills theory as a basis [lS , EH . We derive the relevant Dyson-Schwinger 
equations (or their modification from the pure Yang-Mills theory) for the two-point functions and then consider their 
t-H ' perturbative counterparts. Using techniques based on integration by parts and differential equations, we evaluate the 
OO , noncovariant integrals and present the results for the one-loop, two-point dressing functions. 

The paper is organized as follows. We begin in the next section by considering the first order formalism and deriving 
the relevant field equations of motion. The Feynman rules and general decompositions of the two-point functions are 
obtained in Section 3. In Section 4, the Dyson-Schwinger equations are presented in detail. The one-loop perturbative 
dressing functions are derived in Section 5. The evaluation of the massive noncovariant integrals arising in the one-loop 
■ calculations is presented in Section 6. The results for the two-point functions are presented in Section 7. A summary 
and an outlook are given in Section 8. Various technical details are discussed in the Appendices. 
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FUNCTIONAL FORMALISM 



Throughout this work, we will use the notations and conventions introduced in [14L 1 1 51 ] . We work in Minkowski 
space (until the perturbative integrals are to be explicitly evaluated whereupon we analytically continue to Euclidean 
space), with the metric g^ v = diag(l,— 1). Greek letters (fi,f...) denote Lorentz indices, roman letters denote 
spatial indices and superscripts (a,b,c.) stand for color indices in the adjoint representation. Also, configuration 
space coordinates may be denoted with subscript (x,y,z...) when no confusion arises. The Dirac 7 matrices satisfy 
{7^, 7"} = 2g^ u . The Yang-Mills and quark actions are given by: 
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S q = I d x q(ij Dq + «7 J Dj — m)q 



(2.1) 



(2.2) 
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and the QCD action, Sqcd = Sym + S q , is invariant under a local SU(N C ) gauge transform. In the above quark 
action, q (q) denotes the (conjugate) quark field, the Dirac and color indices in fundamental representation are implicit 
and we have Nf flavors of identical quarks and N c colors. The notation j l refers to the spatial component of the 
Dirac 7 matrices, where the minus sign arising from the metric has been explicitly extracted when appropriate. The 
temporal and spatial components of the covariant derivative (also implicitly in the fundamental color representation) 
are given by: 

D° = d a -igT c a c , 

A = di+igT c Al (2.3) 

where we also rename the temporal component of the gauge field (A 0a ) to a a . The covariant derivative and (in the 
Yang-Mills action) the field strength tensor F are defined in terms of the gauge potential A: 

Kv = d » A i - d »K + 9f abc AlA c v , (2.4) 

where the f abc are the (totally antisymmetric) structure constants of the SU(N C ) group, whose (hermitian) generators 
satisfy [T a ,T b ] = t f abc T c and we use the normalization condition tr(T a T b ) = S ab /2. For later use, the color factor 
(Casimir invariant) for the gap equation will be written Cf = (N% — 1)/2N C . 
Consider now the functional integral: 

Z = J 2?$ exp {iSqcd}, (2.5) 

where T><& denotes the functional integration measure for all fields. Since the action is invariant under gauge trans- 
formations, Z is divergent by virtue of the integration over the gauge group. To overcome this problem we use the 
Faddeev-Popov technique and introduce a gauge-fixing term (for Coulomb gauge: V • A a — 0) along with an associated 
ghost term in standard fashion. At this stage, we also convert to the first order formalism. This entails introducing 
three auxiliary fields: 7?, (f> and t in order to linearize the action with respect to the temporal component of the gauge 
field (a) in local fashion. Classically, if would be the momentum conjugate to A and is here transverse (V • 7? a = 0); 
is a scalar field such that — Vcj> a is the longitudinal component of the conjugate momentum field and r is a trivial 



Lagrange multiplier field. All this is described in detail in Ref. 14] and is not repeated here. Indeed, in this work the 
details of the gauge-fixing and first order formalism are unimportant because the quarks do not directly couple to any 
of these fields (including the ghosts). What is however important later on is that these extra fields will formally enter 
the discussion of the Legendre transform (through partial functional derivatives) which, in principle, gives additional 
terms but which vanish at one-loop order perturbatively. Thus, the reader need only be aware of the existence of such 
fields and be assured that the specific details are not relevant to the present study. 

The full generating functional of the theory is given by the functional integral, Eq. (|2.5|) . in the presence of sources. 
Making the sources relevant to this work explicit and denoting the rest with dots, we have: 

Z[J] = Jv<S>exp^iS QC D + t J d 4 x{x x q x + q x X x +p a <J a + J a -A a + K > a + i? Q -^ + . . .) J. (2.6) 

The field equations of motion (from which the Dyson-Schwinger equations follow) are derived from the observation 
that the integral of a total derivative vanishes, up to possible boundary terms. We use the usual assumption that 
these boundary terms do not contribute [Hj]. For the quark field (we will return to other fields below), we have: 

J V^-^r- exp |i5ym + 1 J d A x [q x (i"/°D 0x + l^Djx - m)q x + x x q x + q x Xx] + ■ ■ • j = (2.7) 

(again the dots represent those source terms that do not play a role here). Using the expression for the components 
of the covariant derivative, Eq. (|2.3p . it follows that 

2?$[(*7°d te + *7 fc V fcx + gT c j°a c x - gT c j k A c kx - m)q x + Xx ] exp {iS} = 0, (2.8) 

where S is the full action plus source terms. 

The generating functional, Z\J\, generates both connected and disconnected Green's functions. However, in practice 
we work with connected two-point and one-particle irreducible n-point Green's functions. The generating functional 
of connected Green's functions is W[J], where Z = e w . We introduce a bracket notation for the functional derivatives 
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of W, such that for a generic source denoted by J a (the index denotes both the type and all other possible attributes 
of the field): 

=<iJ a > ■ (2.9) 

OlJ a 

Converting Eq. (|2.8|) into derivatives of W[J] we obtain : 

(il°d 0x + il k ^kx - m) <iXx> +gT c {^[<ip c x XiXx> + <ip%iXx>] ~ j k [<iJk x ><*Xx> + <iJk x tXx>}} +Xx = 0- 

(2.10) 

We define the generic classical field (we use the same notation for the classical fields and for the quantum fields which 
are integrated over since no confusion will arise) to be: 



1 f 1 sz 

$ a = -J V$»f> aG xp{iS} = - — . (2.11) 



The generating functional of the proper (one-particle irreducible) Green's functions is the effective action, (which 
is a function of the classical fields) and is defined via the Legendre transform of W[J]: 

r[$] = W[J] -i$„J tt . (2.12) 

(We use the common convention that summation over all discrete indices and integration over continuous arguments 
is implicit). This gives: 

sw sr , 

<iJ a >= —— = $ a and <i$ a >= — — = -J a . (2.13) 

OlJ a Ol<P a 

We have used the same bracket notation to denote derivatives of T with respect to fields - no confusion arises since 
we never mix derivatives with respect to sources and fields. Note that care must be taken to observe the correct 
minus signs associated with the Grassmann fields and sources. The equation of motion, Eq. (|2.10|) , in terms of proper 
functions (and from which we will derive the quark gap equation) now reads: 

<iq x >= -i(i-f°dox + il k ^kx ~ vn)iq x + 3 T C 7 ° [a c x q x + <ip c x l X x>] - gT c -f k [A c kx q x + <iJ^ x i Xx >\ ■ (2.14) 

In a similar fashion, the quark contributions to the field equations of motion for <mt x > and <iA$ x > are given by (the 
rest of these equations are simply the Yang-Mills expressions derived previously in Ref. jlij and are not important 
here): 

<«£> = gq x T a 1 Q q x +gTv{T a 1 ° <%%>} + ■■■ , (2-15) 
<iA a ix > = - 5 &TV?x-9Tr {TV <%%>} + •■■ , (2.16) 
where the trace is over Dirac and fundamental color indices. 



3. FEYNMAN RULES AND DECOMPOSITIONS 

Let us now discuss the Feynman rules and decompositions of the Green's functions. The tree-level quark prop- 
agator can be derived directly from the quark equation of motion in terms of connected functions, Eq. (|2.10p . by 
functionally differentiating and neglecting the interaction terms. For the noncovariant case here we obtain (after 
Fourier transforming to momentum space): 

w (o ){k) = _ t fko-j%+m (31) 
— k 2 — m? 

Later on (where appropriate), we will use the usual notation jt= 7°fco — -f l ki. The tree- level gluon propagators needed 
in this work have been derived in [l4T | and are given by: 

fvQ rv n- 
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where tij(k) — Sij — kikj/k 2 is the transverse spatial projector. It is understood that the denominator factors involving 
both temporal and spatial components implicitly carry the Feynman prescription, i.e., 

(3-3) 



(k 2 -k 2 ) (k$-k 2 +l0+) 

such that the analytic continuation to the Euclidean space can be performed. 

The tree-level quark proper two-point function is derived from the quark equation of motion in terms of proper 
functions, Eq. (|2TT4"|) : 

Tf q \k) = l (jk-m). (3.4) 

There are two tree-level quark-gluon vertices, spatial and temporal, again obtained by taking the appropriate functional 
derivatives: 

r& = -gT a i j . (3.5) 

In addition to their tree-level forms, we will also require the general decompositions for the quark two-point functions 
(connected and proper) and the relationship between them. Because we work in a noncovariant setting, the usual 
arguments must be modified to include separately the temporal and spatial components. Starting with the quark 
propagator, we observe that 

W qq (k°,k)~ Jv<$>qqexp{iS} (3.6) 

such that under both time-reversal and parity transforms, the propagator will remain unchanged - the bilinear 
combination of fields is scalar. Since the propagator depends on both k° and k, it has thus four components in 
distinction to the covariant case where there are only two. We thus write 

Wq q (k) = ~=-^ {kol°F t (k) - k a l F s (k) + M(k) + k k i f 1 i F d (k)} (3.7) 

kg — k 2 — TTL 

where all dressing functions are functions of both k 2 and k 2 . At tree-level, Ft = F s = 1, Fd = and M = rn. The 
last term with Fd has no covariant counterpart and in fact will only appear (if at all) at two-loop order and beyond, 
as will be justified below. For the proper two-point function, the same arguments apply and we write 

T qq (k) = i {k 0l °A t (k) ~ k a l A s {k) - B m (k) + k k il °j i A d (k)} (3-8) 

and we will refer to A t ,A s and B m as the temporal, spatial and massive components, respectively. Again the last 
component (Ad) has no covariant counterpart and will only appear at two-loops or beyond. The relationship between 
the connected and proper two-point functions is supplied via the Legendre transform in standard fashion and we have 



r? s (fc)Wft(*) = 1. (3.9) 



In terms of the dressing functions, this gives 

F t = 

F s = 

M = 





(kl-p 


- m 2 )A t 


h 2 A 2 


- k 2 A 2 - 


B 2 m + k 2 k 2 A 2 ' 




(k 2 -k 2 


- m 2 )A s 


h 2 A 2 


- k 2 A 2 - 


B 2 m + k 2 k 2 A 2 ' 




(k 2 -k 2 - 


- m 2 )B m 


h 2 A 2 


- k 2 A 2 - 


B 2 m + k 2 k 2 A 2 ' 




(k 2 -k 2 


- m 2 )A d 


^oA 2 


- k 2 A 2 - 


B 2 m + k 2 k 2 A 2 d 



F d = , rn r„ ' ' ( 3 - 10 ) 



Let us now discuss possible appearance of the genuinely noncovariant term corresponding to the dressing function Ad- 
In this work it will arise, if at all, from the one- loop perturbative form of the quark self-energy (see later for details). 





FIG. 1: Full nonperturbative diagram for the quark self-energy. Filled circles denote dressed propagators and empty circles 
denote dressed vertices. Springs denote connected (propagator) functions, solid lines denote quark propagators and wavy lines 
denote the external legs of the proper functions. 



For now, we can anticipate that since the tree-level quark propagator does not contain a term with fcofci7°7 l and in 
the self-energy loop with two tree-level vertices we only have either two 7 or two 7* matrices together (the gluon 
propagator is cither purely temporal or spatial) , then there is no one-loop contribution that has the overall structure 
7 ° 7 \ This means that Ad = at one- loop order perturbatively. For the rest of the dressing functions, we then get 
the simplified set of relations: 

p (kl-P-m*)A t p = (kl-P-m*)A s M= (fcp 2 - P - m*)B m 

* JL2 A2 _ U2 A2 _ R2 ' S U2 A2 _ U2 A2 _ p>2 ' U2 A2 _ U2 A2 _ R2 \ • ) 

We emphasize that these relations only hold up to one-loop perturbatively — in possible future studies it must be 
recognized that the fourth Dirac structure, 7°7\ may enter in a nontrivial fashion and that the set of relations 
Eq. ([qUO]) should be used. 

Evaluation of the quark contributions to the Waa an d Waa propagators is a direct extension of the results already 
obtained in [l4[ for the Yang-Mills sector. 



4. DERIVATION OF THE DYSON-SCHWINGER EQUATIONS 



We start the derivation of the gap equation by taking the functional derivative of the quark field equation of motion 
(in configuration space), Eq. (|2.14|) . with respect to iq w : 



<iq x iq w >= t(i~f°d 0x +i'y k '\/ kx-m)5(x-w)- / d A yd A z8{x~y)8{x-z) 



r (0)a_<5_ „ r (0)a _S_ ja - 

In the above we have used the configuration space definitions of the tree-level quark-gluon vertices (extracting the 
trivial (^-function dependence in configuration space) and omit those terms which will eventually vanish when the 
sources are set to zero. We use partial differentiation to calculate the terms in the bracket. Notice that simply 
because of the presence of the 7? and <fi fields arising in the first order formalism, we must allow for the additional 
terms so generated (these terms will vanish when we consider the one-loop perturbative case): 

S f 

<ipy%Xz> = - J d 4 vd 4 u <i\ z ix v ><ip a y i(}! a ><%q v iq w %a h u > 

d 4 vd 4 u <iXziXv><'i'PyiKt><tq v iq w i(f>l>, (4.2) 

s c 

<lJjyiX z > = - I d 4 vd 4 u <iXziXv><iJyjiJ b uk><iqviq w iA b uk > 



Siq- 



d 4 vd 4 u <iXziXv><iJyjiK b uk Xiq v iq w ni b uk > . (4.3) 



Inserting the above expressions into Eq. (|4.1[) and Fourier transforming into momentum space we obtain the quark 
Dyson-Schwinger (or gap) equation 

T qq {k) = i ( 7 %- 7 ^ J -m)+ J ^ W r^(A:,- W , W -fc)^Mr| 9CT (c,,-fc,A : -c i ;)IF CT t(fc- W ) 

+ / ^r^ a (A:,- w , w -fc)iF 9 - ? Mr^j w ,-fc,fc- w )^(fc-^) 



FIG. 2: One-loop diagram for the quark contributions to the gluon proper two-point functions. Filled circles denote dressed 
propagators and empty circles denote dressed vertices. Solid lines denote quark propagators and wavy lines denote the external 
legs of the proper functions. 

+ / fcrj;^,^,^ i)^,(a;)r5- TO ( Ul -fc, t - W )< y (fc - (4.4) 

where duj — d A uj / (27r) 4 . The self-energy corrections are presented diagrammatically in Fig. [1] We see that the 
7r and cj> fields do make a contribution thanks to the existence of mixed propagators in the first order formalism. 
But, as emphasized, these contributions will drop out at one- loop order perturbatively because of the absence of 
corresponding tree-level vertices, i.e., that there exist no direct interaction terms in the action between the quark 
fields and the auxiliary fields of the first order formalism. 

We next consider the quark contributions to the proper two-point functions V aa , T a A and Taa- Starting with 
the a equation of motion Eq. (|2. 15[) and following the same procedure as for the gap equation we derive the quark 
contribution to the proper two-point function Y aa in configuration space (trace over Dirac and fundamental color 
indices) : 

<ia a x ia h w > (9 ) = -Tr / d 4 yd 4 zd 4 ud 4: vTf^(z,y,x) <ix y iXu> <Wunviv h w > <iXviXz> ■ (4.5) 



Taking the Fourier transform, we get in momentum space: 



Similarly we 



These loop contributions are shown collectively in Fig. [2] 



(4.6) 

(4.7) 
(4.8) 



r^ (9) (fc) = -'Tr J d U T^{u-k,-u,k)WM^\ q A^^-w,-k)W qq {u-k). 
obtain: 

(i\ _ _ TV / j -p(Q)a/ 7 L\W_ f,V.-n& 



CWfc) - -Try cf W r^ a ( W -fc,-^,fc)^ ? Hr^(a;,fc- W ,-fc)^( W -fc), 



5. ONE-LOOP PERTURB ATIVE TWO-POINT FUNCTIONS 



Let us now consider the one-loop perturbative expansions of the two-point functions derived in the previous section. 
Although so far the formalism has been presented in 4-dimensional Minkowski space, in order to evaluate the resulting 
loop integrals we have to make the analytic continuation to Euclidean space (ko — > iki), where we denote the temporal 
component of the Euclidean 4- momentum k± such that k 2 = k\ + k 2 . Additionally, to regularize the integrals, 
dimensional regularization is employed with the Euclidean space integration measure 



dw A d d w 

* W = (27p+T (5T) 



where d = 3 — 2e is the spatial dimension. To preserve the dimension of the action we must assign a dimension to the 
coupling through the replacement 



2 2 e 

9 -> 9 V > 



(5.2) 
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where fi is the square of a non- vanishing mass scale (which may be later associated with a renormalization scale). 
The perturbative expansion of the two-point dressing function is generically written as: 

r = r(°) + 5 2 r« (5.3) 

where the factor /i £ is included in such that the new coupling and rw are dimensionless. 

In the full nonperturbative gap equation, Eq. (|4.4p . we first insert the various tree-level vertices and propagators 
given by Eqs. (|3.1|) , (|3.2p and (|3.5p . Then we insert the general decomposition of the proper two-point functions, 
Eq. (|3.8p . occurring on the left-hand side of the gap equation, take the Dirac projection, solve the color and tensor 
algebra and lastly perform the Wick rotation. The one-loop temporal, spatial and massive components of the quark 
gap equation in Euclidean space read (recall that Cf = (N 2 — 1)/2N C ): 

Mk) = i- 5 Vc * [ dJ ^ . 2 k ff- 1] A , (5.4) 

kjj \ (lu 2 + m 2 )(k ~ uu) 2 (lu 2 + m 2 )(k - lu) 2 j 
2 en- 1 fj .f 2[k ■ {k - uj)][cj ■ (k - uj)] k ■ lu fc-c3(3-d) 



A 8 {k) = \~g 2 \i e C F — \ duo 



k 2 J 1 (lu 2 + m 2 )(k-Lu) 2 (k-uu) 2 (lu 2 +m 2 )(k-uu) 2 (uu 2 + m 2 )(k - uuf 



(5.5) 

B m (k) =m + rn 9 VC F j dco { {Lj2+ J )Ck _, )2 + ^^(Llu) 2 } ■ (5 - 6) 



•J* u {, 



As mentioned earlier, the possible contribution corresponding to the genuinely noncovariant dressing function does 
not appear at one- loop. To evaluate the integrals occurring in Eq. (|5.5[) . it is helpful to use the identity: 

k-tu= - [k 2 + lu 2 ~ (k - uuf] -k 4 UJ 4 , (5.7) 
which enables us to rewrite A s as a combination of more straightforward integrals: 

(fc 2 + m 2 ) 2 2(k 2 +TO 2 )fc 4 w 4 2efc 2 + 2fc| 



A s (k) = l-sVC f i j dcu J-~- 



oj 2 [(k - uj) 2 + m 2 ]uj 2 uu 2 [(k - lu) 2 + m 2 ]uj 2 [(k — uj) 2 + m 2 ]w 2 
2k-Lu(l-e) 1 m 2 + 3k 2 \ 



[(k - to) 2 + m 2 \u 2 2 [(k — to) 2 + m 2 ]uj 2 



(5.8) 



To determine the quark contributions to the various proper two-point gluon dressing functions given by Eqs. (|4.61 
I4.8[) (presented in Fig. [2]), we again insert the tree-level factors given by Eqs. (|3.1[) , (|3.2p and (|3.5p . solve the color and 
tensor algebra and perform a Wick rotation. The one-loop integral expressions are: 



k 2 T yL > Jk) = ffN f 2 / dui — - — - — - — ^ttt, — " \, , (5.9) 
" "(*) v ' 1 ] (lu 2 + m 2 )[(k - uj) 2 + to 2 ] ' v ; 

wS^M - *» ,^ + ,:^;^ + mT (»■«>) 

(5.11) 

where , and T *AA(q) are ^ ne transversal and longitudinal components of the proper two-point function r^J^ 

(given by Eq. (14. 8p ). respectively (see Ref. [l4[ for details of decomposition). 



6. NONCOVARIANT MASSIVE LOOP INTEGRALS 



In the one-loop expansions from the previous section, there are two types of integrals arising: those which can 
be solved using standard techniques (such as Schwinger parametrization, Mellin representation - for details, see 
Appendix [A| , and those which require a more esoteric approach. In this section we study the latter variety using a 
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technique based on differential equations and integration by parts developed previously [1 5I | . We will consider the two 
integrals: 

A m (kl,k 2 ) = ( — — 21^2' (G-l) 

J uj z [(k — uj) z + m z \uj z 

Amik^jk 2 ) = / 2777 yT~ 1 21 -i ■ (6-2) 
J uj z [(k — uj) 2 + m z \uj z 



A. Derivation of the differential equations 



Let us first write Eqs. (|6.ip and (|6.2p in the general form (n = 0, 1) 

duj uj™ 



uj 2 [(k — to) 2 + m 2 ]uj 2 



(6.3) 



In this derivation k\ and k 2 are treated as variables whereas the mass, to, is treated as a parameter. The two first 
derivatives are: 



k 

4 dki 

gjn 

ok k 



uj 2 [(k — uj) 2 + m 2 ]uj 2 



-2 



(fc — Uj) 2 + TO 2 

fc-(ife-cZf) 

u> 2 [(k — uj) 2 + m 2 ]uj 2 I (k — uj) 2 + m 2 



dbJ UJ™ 



There are also two integration by parts identities: 



: 




a 



n+1 



dbJ 



duJi uj 2 [(k — uj) 2 + m 2 ]uj 2 



du>i uj 2 [(k — uj) 2 + m 2 ]uj 2 
Adding these two expressions gives 

doj uj™ 



uj 2 [(k — uj) 2 + m 2 ]uj 2 

duj uj™ 
uj 2 [(k — uj) 2 + m 2 ]uj 2 



n + 1 - 2-t- - 2 



d-2-2- 



uJiiuJi — k^) 
(k — uj) 2 + TO 2 

uj - (uj — k) 
(k — uj) 2 + TO 2 



= 



uj 2 [(k — uj) 2 + m 2 ]uj 2 
Expanding the numerator factor, we can rewrite Eq. (|6.8p as 

du> u>2 



d+n-3-2 



uj ■ (uj — k) 
(k — uj) 2 + TO 2 



' uj 2 [(k — uj) 2 + m 2 ]uj 2 
Combining this with Eq. (|6.6p then yields: 
duj 



fc 4 (fc 4 - UJj) 

uj 2 [(k — uj) 2 + m 2 ]uj 2 \ (k — uj) 2 + rri 



duj 



k- — uj z + TO 
(k — uj) 2 + TO 2 



(6.4) 
(6.5) 

(6.6) 
(6.7) 

(6.8) 
(6.9) 



duj uj\ 



k 2 + to 2 J [(k — uj) 2 + ' 



duj uj"2 



uj 4 [(k - uj) 2 + to 



uj 2 [(k — uj) 2 + m 2 ]uj 2 
2 



2k 2 



k 2 



duo UJ™ 



:(n + d-4)-n+l 



uj 2 [(k — uj) 2 + TO 



212 ' 



(6.10) 



This leads to the temporal differential equations for A m and 



&4 



dA m 

<%4 



1 + 2 



(d-4)k 



k 2 



Ar. 



duj 



k 2 + to 2 



k 2 + to 2 J [(k — uj) 2 + 



duj 



uj 4 [(k - uj) 2 + m 2 



duj 



uj 2 [(k — uj) 2 + TO 2 ] 2 ' 



dk,A 



M-3)k 2 
k 2 + to 2 



A, 



d uj UJ4. 



k 2 + to 2 J [(k — uj) 2 + m 2 ] 2 uj 2 



duj UJ4 



UJ 4 [(k - uj) 2 + TO 2 ] 



duo uj^ 



uj 2 [(k — uj) 2 + TO 2 



(6.11) 



(6.12) 
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In the same manner, we derive the differential equations involving the spatial components: 



kj 



dA„ 



2-d + 2 



(d-A)k 2 
k 2 + m 2 

tibJ 



A,, 



2k 2 



k 2 



[(k — lu) 2 + m 2 ] 2 LU 2 



LU 4 [(k - Lu) 2 + TO 2 



Lu 2 [(k - lu) 2 + m 2 ] 2 
(d - 3)k 2 



(6.13) 



2-d + 2 



-4 



2k 2 



dbJ UJ4 



k 2 + m 2 J [(k — lu) 2 + m 2 ] 2 LU 2 



duj LU4 



Lu 4 [(k - lu) 2 + m 2 



dbJ LU4 



uj 2 [(k — lu) 2 + m 2 



(6.14) 



It is in fact possible to write down a mass differential equation, but in the light of the method presented here this 
would not bring any new information. However, this third derivative is important because it provides an useful check 
of our solutions (for a detailed discussion, see Appendix [B| . 

At this point, let us discuss how the differential equations for the massless integrals considered in Ref. [l5[ are 
regained. In the massless limit, there are potential ambiguities arising in the integrals appearing in Eqs. (|6.11l - l6TT4"]) . 
because in part, the limits m —> and e — ► do not interchange. Let us start by considering the following integral 
given by Eq. (|A.16[) (similar arguments apply to all the integrals appearing in the differential equations): 



I = 



[m 2 ]- 1 - £ r(i-e)r(l + e) 



F\ l,l + e;3/2-e; 5 



(6.15) 



to 2 [(k- lu) 2 + m 2 ] 2 (4tt) 2 - £ r(3/2 - e) 
It is useful here to invert the argument of the hypergeometric with the help of the formula (see, for instance, Ref. [19|): 

r(c)r(&- a) 



2 F 1 (a,b; c; t) 



T(b)T(c-a) 
T{c)T{a-b) 



{-t)~ a 2 Fi [a, 1 - c + a; 1 - b + a: 



1 



{-t)-° 2 F 1 [b,l-c+b;l-a + b;- 



Then we have: 



I = 



1 [m 2 ]- £ T{e) 



F^l^ + s-A-e;-^ 



T(a)T(c - b) 

[k 2 ]- 1 - £ T{\~e)T(-e)T{l + s 



(6.16) 



(4tt 



,2-e 



r (2 - 2 ^) 



1 m 2 
■2F1 l+e, - + 2e;l + e;- — 

2 k 2 



(6.17) 



In the expression above, the problem of the non-interchangeable limits is seen explicitly in the first term. However, 
when all the integrals occurring in the various differential equations are put together, such terms explicitly cancel and 
only the second term of Eq. (|6. 1T[) (which leads to the correct massless limit) contributes. 

Returning to the differential equations, we evaluate the standard integrals in terms of e (see Appendix [AJ and with 
the notation x — k 2 , y — k 2 we obtain for A m : 

OAm 



2x- 



dx 



x + y + m 2 



2y- 



dy 



l-2(l + 2e) 

r(- g )r(i + e) 

r(2-e) 
-l + 2£-2(l + 2e) 

r(- e )r(i + £ ) 



,21-1- 



A r , 



2F1 ( 2,l + e;2-e 

y 



x + y + in 1 



(4tt) 2 - e 
x + y 

A — 9 ^ 



y + m 
x + y + m 2 

Y 2 F 1 (l,l + e;2-e 

y 



X 2 Fx fl,l + e;3/2-e;— \ 
V m 2 

x + y 



(6.18) 



(4tt; 



2-e 



x + y + m? 



Y{2-e) 
and for the integral A 4 = k^A T 
dA m 



Fx 2, l + e; 2-e; 



x + y 



Y 2 F X 1, l + e; 2-e 



X 2 F X (l,l+e;3/2-e;--^ 
x + y 



(6.19) 



2a;- 



dx 



■l-4e- 



x + y + m 2 



A, 



[m 
"(4^ 



21-1-6 



)2-e 



y + m 
x + y + m 2 



X 2 F X (l,l + e;3/2-e; 



m - 



10 



Y 



2y- 



dy 



2-e 
-l + 2e-4e 
Y 



Fi 2,l + e;3-£;- 



x + y 



y 



x + y + 
2 F l ( 2, 1 + e; 3 - e 



A m — 2 
x + y 



in 



{l-e) 2 Fi [ l,l + e;3-e; 

21-l-e 



(4tt) 2 6 \x + y + m 2 
-{l-e) 2 F 1 ( l,l + e;3-e; 



Xa-Fi (l,l + e;3/2-e;— 2j 



a; + y 



where 



Y _ r(i/2- £ )r(i + E ) y= r(i-e)r(i + e ) 



r(3/2-e) 



r(2- £ ) 



(6.20) 



(6.21) 



(6.22) 



B. Solving the differential equations 

Let us first consider the integral A m . By the same method as in [l5j], we make the following ansatz: 

A m (x,y) = FA m {x,y)GA m (x,y) 



such that 



2x 

FAmZx 
FAmty- 



dF Ar , 
dx 

dF A „ 
dy 

dx 



1- (2 + 4e)- 



x + y + m 2 
-l + 2e- (2 + 4e) 



Fa* 



l) 



= 2 



,2i-l-e 



(47r) 2 e |_ x + y + m 

r(- e )r(i + e ) 



x + y + m- 

2 



^±^ M (l,l + £ ;3/2- £ ;-4 



9y 



r(2- E ) 

„2l-l-e 



= 2 



2 Fi ( 2,l + e;2-£ 



a; + y 



Y 2 F 1 l,l + e;2-e; 



+ 



(47r) 2 6 \ x + y + to 2 

r(- E )r(i + e ) 



X 2 Fi (i,i + e ; 3/2 -e;--% 



F x 2,l + e;2- £ ;- 



x + y 



Y 2 Fi l,l + e;2-e; 



£ + y 



x + y 



(6.23) 

(6.24) 
(6.25) 

(6.26) 

(6.27) 



r(2 — e) \ y \ 

By inspection, it is simple to determine the solution for the two homogeneous equations, Eq. (|6.24p and Eq. (|6.25() : 

F Am (x, y) = x^y-V^ix + y + m 2 )- 1 " 2 *. (6.28) 



Since the mass m is treated as a parameter, the (dimcnsionful) solution, Eq. (|6.28p . may have an integration constant 
proportional to [m 2 ] _1_£ . However, returning to the original equations, Eq. (|6.24[) and Eq. (|6.25[) . we see that the 
only consistent solution is the one for which this constant vanishes. 

For the function Gatu we make the following ansatz, which will be verified below (z — x/y): 



G Am {x,y) = G° Am (x,y) + G Am (z). 



(6.29) 



The component G° Am (x, y) can be found by adding the differential equations Eq. (|6.26p and Eq. (|6.27[) . which lead to: 



9G° Am dG A 

x _Arn + y A 



dx 



In 



1 + ^ + 1 



dy (4tt) 2 5 y y/x(y + m 2 ) ^ ^l+ ,J f 



- 1 



0{e) 



(6.30) 



Because the function G Am is multiplied by the function F Am (which does not have an e pole), the term of order 0(e) 
will not contribute. The solution of this equation is: 



G°Am( x i y) 



(4tt 



.2-6 



1 + — In 



- Iny 



0(e) } +d. 



(6.31) 
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Before we proceed to determine GAm, we justify the ansatz for GAm, given by Eq. (|6.29|) . First we observe that: 



2xF A r 



9Gat, 



= 2zF A , 



9Gat, 



+ 2xF Ar 



dG° A , 



OX OZ ox 

Then we subtract the above equation, Eq. (|6.32[) , from Eq. (|6.26|) . This gives: 



(6.32) 



dz 



1 [m 2 Y 



F Am (4tt) 2 -^ 

r(- £ )r(i + £ ) 

T(2-e) 



y + to 

X + y + 777 2 



X 2 Fi (l,l + e;3/2-e;- 



Y 2 F 1 l,l + e;2-e: 



x + y 



F! 2, 1 + e; 2 - e; 



dx 



(6.33) 



Evaluation to the first order in e is straightforward and we see that the right hand side of the above expression is only 
a function of the variable z. This allows us to write down a differential equation for GAm(z) in the form: 



dG Ar , 



1 



1 I 1 



(47T) 



2-e 



/Z £ 



- -7 + hiTO 2 +C3(e) \ , 



from which we get immediately 



GAm(z) = - 



1 I 1 



(4tt 



(2-e 



<Z I £ 



- -7 + lnm^ +0(e) } + C 2 



(6.34) 



(6.35) 



Returning to the original differential equations (|6.24I - HT2"T)) with the function G(x, y) = G Am (x, y) + GA m ( z ), we see 
that the only consistent solution is the one for which the overall constant C\ + C 2 vanishes. 

We may now put together the solutions Eqs. (|6.28[) , (|6.31[) and (|6.35p and write for the function A m : 



A m (x,y) 



(x + y + 771 ) 



2\-l-e 



(4tt 



,2-e 



h 27 + 2 In 




We see that for to 2 — we regain the result from [l5| and that the singularities are located at x + y + m 2 = 
(with to 2 , y > 0). Two more useful checks arise from the study of the power expansion around x = and the mass 
differential equation (for details, see Appendix |B|) . 

We now proceed in the same fashion to determine the function A m (x,y) — F-r {x, y)G-^ m (x, y). The resulting 
partial differential equations are in this case: 



2.r 

2y 



FAm^ 



FArJV 



dF-j 

Am 




-1 - 


dx 




dF-j 

Am 

dy 




-1 + 


dGAm 


[to 2 ] 


dx 


(4tt)< 






Y 






2-e 


dGAm 


Jm 2 }- 


dy 


(4tt) 5 






Y 




+ 


2-e 



x + y + to 2 

. a. y 



X + y + 777 

y + m 2 



^Am ' 



x + y + 777 2 

2 .F L (2.1 • : :', 

y 

x + y + to 2 



X 2 F 1 (l,l + e;3/2-e; — \ 

\ 771 

x + y 



(l-e) 2 Fi l,l + e;3-e;- 



X 2 F 1 1,1 + e; 3/2 - e; -- ^ 



Fi 2,l + e;3-£; 



x + y 



(l-e)aJi l,l + e;3-e; 



x + y 



to 2 



with X, Y" defined previously. The solution to the first pair is 

2 v 



F-jJx, y) = x-^y-^ix + y + to 2 )- 2 -. 



(6.37) 
(6.38) 

(6.39) 

(6.40) 
(6.41) 



For brevity, in the above expression and also in the derivation of the function G^ m — G^ + G^ rn (the analogue of 
GAm) we omit the constants of integration - they vanish as in the case of the functions FAm and GAm- 
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As before, for G-j (x, y) we make the ansatz: 



G Am( x > y) = G L.( X ' V) + ^ImW' 



Am 

In the limit e — > 0, the component (x, j/) is determined from the differential equation: 



(6.42) 



dG^L dG^L 

~ Am , Am 
a " a 

ax ay 
The solution of this equation is: 



(4tt) 2 



+ y + m 2 ^Jy + 



;ln 



1 + ^ + 1 



1 + J21-1 



0(e) 



(6.43) 



-iLi 2 



(4tt 



|2-e 



iln 



In 



l + ^ + i ^. 



zVI+ 1 
1 



1 - iyz 
1 + iJz 



1 + s^-iJz 
v v 



iLi 2 



1 + ^+z^i, 



1 + zyl 

1 - iJz 



1 + - t^fz 



i + ^ + ^Vi, 



where hi2(z) is the dilogarithmic function [2C 



Li 2 (z) 



ln(l - t) 



dt. 



0{e) 



(6.44) 



(6.45) 



/o t 

As before, we check that the ansatz for G^ m (a;, y) given in Eq. (|6.42[) is correct and derive the differential equation 
for the function G-^ m , in the limit e — > 0: 



1 



(47r) 2 - £ [z+1 
The result we leave for the moment in the form: 
1 



[ln(l + z) - In z - 2 In 2] + 0(e) 



G~ Am {z) 



(4tt) 5 



-4 In 2 arctan(v / i) 



ln(l +t) 



dt 



o V£(l+t) 



hit + 0(e) 



(6.46) 



(6.47) 



With the solutions, Eqs. (|6llj) . (jQ4|) and fOT]) . after some further manipulation we can write down the following 
simplified expression for the integral A^: 



A 4 m ( X ,y) = fc 4 V " - ? ^ V ' ' "' " ' <j 2 b (— ) arrt.an + 91t, | V . V larctau 



( 47r )2-e 

dt 



1 + 



y 



with the integral 
dt 



In 1 + 1 



0(e) 



(6.48) 



o Vi(l+t) 
\-i-Jz\f 1 



In I i + =7 rln2-ilnf ^— ^£ ) 2h:2 



■ In 



, 2 

y 

iJz 



l + lx/Z 



In + + 1 



■ In 



1 + 



1 



W2 



— z In (y« — «) 
+ iln{Jz + i) 



In 2 + In (1 + z H ) - In (1 - iJz) In (Jz - i) 

V 2 

In 2 + In (1 + z H ) - In (1 + i-y/i) ~ ln (V£ + «) 



1 



+>.(5-|V5)+»Li2(5 + 2^ 



iLi 2 ■ — 7= + «Li 2 



iLi 2 



M 1+^ + ^ 



zLi 2 



, \- r -f+Jz 

i + Jz 



-,w'l + ^ + 
-« + \/i 



-i + sJZ 

-j~z 



-l+y/Z 

I + Jz 



iLi 2 



i + Jz 



(6.49) 
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We see that for to 2 = we get the correct limit for the function A^. We also mention that the singularities are 
located at x + y + m 2 = and the apparent singularities at z = — 1 (i.e., x + y = 0) in the expression Eq. (|6.48[) are 
canceling out. This can be easily seen by making a series expansion of Eq. (|6.48|) around z = — 1: 



At 



1 



z + l 



m - 



6i) +0((z + lf 




Again, the result Eq. (|6.48p has been checked by performing an expansion around x 
differential equation (see Appendix |B|) . 



(6.50) 



and by studying the mass 



7. PERTURB ATIVE RESULTS IN THE LIMIT e -> 



We can now collect together the results and write down the one-loop perturbative expressions for the two-point 
functions. In Eqs. (|5.4[) . (|5.6[) and (|5.8[) we insert the corresponding integrals (derived in the previous section and 
Appendix[XJ and find for the temporal, spatial and massive components of the quark gap equation, in the limit e — > 0: 



A t (k) = 1 






C F9 2 I 
(4^) 2 " e 1 




2 

TO 

— m 


A s (k) = 1 






<W» I 
(4^) 2 - £ 1 




2 

TO 

— m 






/ m 2 


- 4— + 

V fc 2 





/fc 4 



l)ln(l + 4)+^)^ (7-1) 



9 

771 z 





B m {k) 

_Ci 

' (47r) 2 - e |e M *"V ~ A: 2 \ , ,'| . £ ._ [ / " V k' 1 J "' V nr 

where the function f m (x,y) is given by (z — x/y = k 2 /k 2 

2 ( m 2 \ 2 /V 1 + ^ _ + 1 \ 
f m (x, y) = —= In — ) arctan (\fz) H p= In — arctan 



/o Vt(l + rf) 
The last integral has been rewritten using the identity: 





(7.2) 



(7.3) 



ln( i + ^+^l. (7.4) 



dt , / to 2 \ r df , / m 2 \ . , 

- / -p; 7 In = / -f; r In 1 + art + — . 7.5 

zj Vt(i + t) V W \/t(l + zt) V 2// 

As a useful check, we can set m = and show that the results for the temporal and spatial components are in 
agreement with the calculation performed independently using the one- loop massless integrals derived in Ref. [l5| • As 
has been shown in the previous section, in the noncovariant integrals the singularities appear at x + y + to 2 = 0. It is 
easy to see that the standard integrals also have the same singularity structure. Because of the absence of singularities 
in both the Euclidean and spacelike Minkowski regions, we can see that the validity of the Wick rotation is justified. 

Having calculated the dressing functions for the quark proper two-point Green's function, we are now able to discuss 
the structure of the propagator. In Eq. (|3.1ip we first analyze the denominator factor. Let us denote (in Euclidean 
space): 

D(k) = klA 2 (k) + k 2 A 2 s (k) + B 2 m (k). (7.6) 
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Inserting the expressions from Eqs. (|7.ip . (|T. 2[) and (|7.3|) into the above equation, we have: 



D(k) = k 2 + m 2 \l + 6 



g 2 C f 

2-e 



(47T) 



In- 



a 2C F <? 
m ) 



, n 2 2\ 2Cf3 
+(fc + m ) 



(47T) 



2-e 



fc 2 



3 - — In 1 + - 4W 1 + — In 



We define the renormalized mass, trr, via: 



2 = Z^m 2 ^ with Z 2 „ = 1 — 6 



g 2 c F 

( 47r )2~ £ 



fc 2 




.(7.7) 



hi 



m 



(7.: 



The expression for D(k), Eq. (|7.7[) . then contains explicitly the overall factor k 2 + m 2 R: meaning that the simple 
pole mass of the quark emerges, just as it does in covariant gauges. The singularity structure of the remaining part 
is such that there are no non-analytic structures for spacelike or Euclidean momenta. Moreover, we see that the 
renormalization factor, Z mi which should be a gauge invariant quantity (it defines the physical perturbative pole 
mass) agrees with the result obtained in covariant gauges [2l| . 

Because of the Dirac structure, it is more convenient to write the quark propagator in Minkowski space. We have 
shown that the analytic continuation of the functions A t ,A s , B m back into the Minkowski space is allowed and this 
enables us simply to write (note that also in D(k) we must also analytically continue k\ — > — feg): 

Wg q (k) = i {fk a A t {k) - fk t A s (k) + B m (k)} D-\k). (7.9) 

Inserting the denominator factor, Eq. (|7.7p . in the limit e — » and replacing the mass with its renormalized counter- 
part, the above expression gives: 



W- m {k) 



(jt+m R ) 



1-Cf 



9 



(47T) 



2-e 



finite terms 



We can thus write down for the quark propagator: 



W- qq {k) = - 



i(jt + m R ) 



fcp — k 2 — m 



finite terms 



and identify the renormalization constant (omitting the prescription dependent constants) 

9 2 C F 



Z 2 = l- 



(4tt 



,2-e 



- - 7 



(7.10) 



(7.11) 



(7.12) 



Turning to the quark loop contributions to the gluon two-point proper functions, in evaluating the integral structure 
of Eqs. (|5.9[) . (|5.10[) and (|5.11[) we observe the following relations (in Euclidean space): 



(7.13) 



where the integral I(fc|, k 2 ) reads (using the results of Appendix VK\ , as e 
N f 



I(kik') 



+ 



( 47r )2-e 

2 



in 



2 

_ 3 

In- 



In- 



0(e) 



k 2 



10 



4m 2 




(7.14) 



The above integral agrees with the results obtained in covariant gauges (see for instance |21[ ) . This is hardly surprising, 
since at one-loop level the quark loop as a whole is unchanged from its covariant counterpart — what is different is that 
the various degrees of freedom (temporal and spatial) are being separated into the corresponding proper two-point 
functions, i.e., Ta<j and r^. 
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The one-loop gluon propagator dressing functions we construct by writing D = + g 2 D^K As mentioned 
previously, in the first order formalism we have to account for the presence of the additional 7?, (j) and ghost fields and 
the corresponding propagators (for example Da-k)- Whilst the quarks only contribute to three of the gluon proper 
two-point functions (r^, T Aa and T aa ) at one-loop, there will be contributions to many more of the various connected 
(propagator) two-point functions. The relationship between the connected and proper gluon two-point functions in 
the first order formalism is detailed in Ref . [l5| • The full set of quark contributions to these gluonic type propagators 
is: 



2 



D 



% q) (k) = = D^ {q) {k) = -D^ {q) (k) = I(km (7-15) 



At this point we are able to identify the first coefficient of the perturbative /3-function. As is well known in 
Landau gauge, a renormalization group invariant running coupling can be defined through the following perturbative 
combination of gluon and ghost propagator dressing functions [22j: 



9 2 D A aD 2 c ~ g 



1 , 9 2 1 (UNc 2N f \ 
16tt 2 e \ 3 3 ) 



(7.16) 



At one-loop in perturbation theory, the coefficient of the 1 /e pole above is simply minus the first coefficient of the 
/3-function (/3q = — llN c /3 + 2iV//3). By inspecting the relations Eq. (|7. 15[) and those obtained in Ref. [l5| for the 
propagators Daa and D c , we see that the same result is achieved in Coulomb gauge. In Coulomb gauge, a second 
renormalization group invariant combination of propagators appears and is given by g 2 D aa Again, combining our 
results Eq. (|7. 15|) and those obtained in [15j, we see that indeed the coefficient of 1/e agrees with this. 



8. SUMMARY AND OUTLOOK 



In this paper, the quark contributions to the Dyson-Schwinger equations of QCD have been derived within the 
Coulomb gauge first order formalism and perturbative results have been presented. The set of Feynman rules has 
been derived and the general form of the two-point functions have been established. The quark gap equation and the 
quark loop contributions to the Dyson-Schwinger equations concerning the gluon proper two-point functions have 
been explicitly derived. A one-loop perturbative calculation has been performed, for the quark gap equation, as well 
as for the quark contributions to the gluon proper two-point functions. The required noncovariant massive integrals 
have been obtained, using techniques based on differential equations and integration by parts. The various two-point 
dressing functions and propagators have been evaluated in the limit e — > 0. The validity of the analytic continuation 
between Minkowski and Euclidean space has been verified. The quark mass and propagator have been renormalized 
and it is seen that the one-loop result for the gauge invariant quark mass renormalization coefficient agrees explicitly 
with the result obtained in linear covariant gauges. The correct one-loop coefficient for the [3 function has also been 
obtained. 

The natural continuation of this work is the perturbative evaluation of vertex functions of the theory. The Mellin- 
Barnes parametrization or perhaps a generalization of the differential equation method to the three-point integrals 
would be possible ways to proceed. Also, the construction of scattering matrix elements would be another interesting 
topic. 
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APPENDIX A: STANDARD MASSIVE INTEGRALS 



Consider the integral: 



J m {k 



,21 M 



[(k 



(A.l) 
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In the case /u = v = 1 this gives the scalar integral associated with, for example, the fermion loop in quantum 
electrodynamics 1211 . We present here a method to evaluate such integrals for arbitrary denominator powers (developed 
originally in Ref. [23]) and generalize to the various additional noncovariant integrals. 

We start by writing the Taylor expansion of the massive propagator in terms of a hypergeometric function in the 
following way: 

1 1 / m 2 \ 

[ui z + m z \i J - [w 2 J^ \ uj z J 

Now, the idea is to use the Mellin-Barnes representation of the hypergeometric function iF (^l; z): 

zoo 

1 F (fi;z) = f ^-^- J ds{- Z yT{-s)T{n + s), (A.3) 

—zoo 

where the contour in the complex plane separates the left poles of the T functions from the right poles. A first advantage 
of this representation is that the "mass term" gets separated from the massless propagator and the remaining integrals 
can be calculated with the Cauchy residue theorem, as we shall see below. Also, we note that the results can be written 
as a function of either k 2 /m 2 , or m 2 /k 2 (expansions thereof are of interest in studying various momentum regimes). 
This we do by using the well-known formulas of analytic continuation of the hypergeometric function (for an extended 
discussion, see [Hj]). 

Applying Eq. (|A.3[) to the massive propagator we can rewrite the integral J m (k 2 ) as: 

zoo 



Inserting the general result for the massless integral (an explicit derivation can be found in Ref. [Ta]): 

duj _ [k 2 f-^- v - e T(fi + v + e - 2) T(2 - fi - e)T(2 -v-e) 



[w 2 H(fc- uj) 2 Y (4tt) 2 - £ T{h)T{u) r(4 - fi - v - 2e) 

we get for the integral J m : 



(A.5) 



L __J / / J„JJ. / ! 

(47r) 2 " e (2m) 2 T(ji)T{u) 



J^k 2 ) = 1 > e <0 II dsdt(-) T(-s)T(-t)T(2- S -fi-sn2- S -u-t) 



T(fi + V + s + t-2 + e) 



r(4 - 2e - fi - v - s - 1) 



With the change of variable t = 2 — e~ fi — v — u — s (for such a replacement, the left and right poles of the V function 
are simply interchanged and therefore the condition of separating the poles is not contradicted) we obtain: 



v.oe 



(4tt) 2 - £ (2m) T(fi)T(v) J \k 2 J T(2 - e + u) 

— ZOO 

too 

^J— [ dsT(-s)T(2-£- fi- s)T(-2 + e + u + fi + u + s)T(fi + u + s). (A.7) 
m) J 



(2m) 

— zoo 

To evaluate the integral over s we use the Barnes Lemma: 



1 f j n i nw, m,/ s^tj ^ T(a + c)T(a + d)T(b + c)T(b + d) 

1 ds T(a + s)T(b + s)T(c - s)T(d - s) = — '—± -(— ^ (A. 



(2m) J T(a + b + c + d) 

—too 

and for the integral Eq. (IA.7|) it follows immediately that 

\ m 2i2-n-v-e x 1 > / m 2\ « r(-u)T(fi + u)T(v + u)T(p + v - 2 + £ + u) 



Jm{k ] - (4tt) 2 ^ (2m) T(n)T(u) J dU U 2 J T(n + v + 2u) 
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Closing the integration contour on the right we have: 

(47r)^ 6 (2m) T{n)T[y) f~^\ k J 3- r(/i + z/ + 2j) 



With the help of the duplication formula 



r(2z) = 2 22 - 1 vr- 1 / 2 r(z)r [* + -), (A.n) 



we can rewrite J rn as: 

[m 2 ] 2 -^-"- 6 + i/ - 2 + e) 



Jm{k ) 



(4tt) 2 - £ r(/x + 1/) 

V f-^-V 1 i rfc + jW + ^rfc + ^ + e + j) r($) r(*±p) 
^ ' fc 2 y 2 2 Jji r( M ) r(i/) r( M + ^-2 + £ ) r(^ + j) r(^±^ + j) ' 1 ' 1 



The sum is clearly a representation of the hypergeometric ^F 2 {a^ 6, c; d, e; z) (see (19() and we finally obtain: 

J ^ k ) - d) 2 - g 3F2 ^-^ + - 2+£; V'^^ ; -4^j- (A - l3) 

A trivial computation shows that the result Eq. (|A.13|I is consistent with the known results in the limit m = 0. All 
we have to do is to invert the argument of the hypergeometric according to the formula (see, for example, Ref. |24|) 

Fin n n -h h ■ A - T ^ T ^ f r ( Q l) F ( Q 2 " gl)T(g3 ~ gl) / a -a, 

3 ^ 2 (ai,a2,a3,0i,02,zj - — — — : — — - — - < -— — — (-z) 

r{ai)T(a 2 )T(a a ) { T{bi - ai)T(b 2 - ai) 

X3F2 ( ai,Oi - 61 + 1, a! - b 2 + 1; ai - 02 + 1, ffli - 03 + 1; - 
V z 

r(a 2 )r(ai - a 2 )r(o 3 - a 2 ) , ._„, „ / . , 1 t , 1 , 1 , -, 1 

-(-2) a2 3 F 2 a 2 ,a 2 -61 + 1,02-62 + 1; -01 +a 2 + 1,02 -03 + 1; - 



r(6i-a 2 )r(6 2 -a 2 ) v y - - - ^ - — > z 

r(a 3 )r(ai - a 3 )r(a 2 - "^-a^ ? ^_ hi + ^ ^ _ ^ + 1; ^ + ^ + ^ _ fl2 + fl3 + r 1 



r(6i - a 3 )r(6 2 - a 3 ) \ " " " " z 

(A.14) 

For integrals with different type of denominator factors, similar calculations bring us to the following results: 
duj _ [to 2 ] 2 -^- £ T(2 - fi - e)r{ji + is + e - 2) f k 2 

[io 2 Y[{k - uj) 2 + m 2 Y ~ {att) 2 -z r(i/)r(2 - e) 2 i y 1 ^ i + " + £ - ; ~ £; ~;^2 

(A.15) 

tfw [m 2 ] 2 -^-^- £ r(|-/i-e)r( A i + z/ + £-2) / j?\ 



/" ^ _ [m 2 ] 2 -^- g r(|- M -e)r( M + ^ + g -2) / _ 9 o/o_ 

7 [£ 2 M(fc- w ) 2 + TO 2 ]^ (4^) 2 - I»r(3/2 - e) 21 ^^"^ e; 



9 

777 z 



(A.16) 



This method can also be applied to integrals with more complicated numerator structure. 

For completeness, we also show the first order e expansion of the integrals arising into the one-loop perturbative 
expressions considered in this work: 



(c 2 + to 2 ) [(k c) 2 + to 2 ] ( 47 r) 2 - \~e ~ ^ + 2 ~ V 1 + J^ 1 * [ TT^Z[ ' + ^ ' " ' ' ' ' ? ' 



A' 2 



fm 2 l- £ f 1 / to 2 \ / fc 2 . , . 

7 + 2- 1 + tt In I + -5 +0(e) , (A.18) 



uj 2 [(k - uj) 2 + m 2 ] (47r) 2 ~ e (.£ V ^ 2 







, [to 2 ]^ £ (l/l \ 1 lm 2 1/ to 2 \ 2 , / k 2 , , 



w 2 [(fc-w) 2 +TO 2 ] 1 (4tt) 2 ~ £ 2 Ve 7 2 fc 2 2 V k 2 



771 z 
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(A.19) 

/ , 2Uh % 21 = ^h-2, + 8 -2,I^J^ +1 ) + Q {s) ). (A.20) 
J uj 2 [(k - lu) 2 + m 2 } (47r) 2 - £ | £ y fc 2 ^ / 1 + j S ?_ 1 J W j 



APPENDIX B: CHECKING THE NONSTANDARD INTEGRALS 



One way to check analytically the results for the integrals A m and A^, Eq. (|6.36p and Eq. (|6.48p . respectively is to 
make an expansion around x — and evaluate the resulting integrals with the help of the Schwinger parametrization. 
Let us consider then the integral A m , originally defined in Eq. (|6.1|) . Using Schwinger parameters [25l j , we can rewrite 
the denominator factors as exponential functions to give: 



A, 



dadftd'y / dus exp 



I -(a + ft)u\ + 2/3/c 4 w4 -0kl-{a + P + ~/)lS 2 + 2(3k-uj - (3k 2 - 0m 2 }. (B.l) 



Applying similar reasoning as in Ref. [la ] , we come to the following parametric form of the integral (recall that . 



An 



(x + y + m 2 )" 



(4tt; 



2-6 



-r(l + e) / d/3 



1-0 



da 



(a + /3)V2 



a(3 



(3(1 - (3)y + (3m 2 



(a + (3) (x + y + to 2 ) 



x + y + to 



-l-e 



(B.2) 

For general values of x, the integral above cannot be solved because of the highly nontrivial denominator factor. Since 
there can be no singularities at x = (this would invalidate the Wick rotation which, as discussed in the text, does 
hold here), we make an expansion around this point and then integrate. To first order in powers of x we have: 

-l-e 



Ar. 



*o (x + y + m 2 y 



(4tt) 2 - £ 
to 2 + y(l - (3) 



1-/3 



r(l + e) / dp / da(a + (3)- 1/2 



(3 



Til + y(l - (3) 



ft 



-2-e 



m z + y 

After performing the integration we get: 

, , v x^O (x + y + m 2 )- 1 - 6 
A m (x,y) = 



aft 



(m 2 + y)(a + (3) 



ft 



to 2 + y(l - (3) 



(to 2 + y)" 2 



Tn + y 



(1 + e)x + 0(x 2 ) 



(B.3) 




r(l+e) 2 F 1 _ e ,2 + e;l-e; 



to^ + y 



+ 0(x 2 )+0(e) 



(B.4) 



In the above formula, we isolated the hypergeometric term and evaluate the e expansion separately. In order to do 
this, we have to differentiate the hypergeometric function with respect to the parameters. In general, differentiation 
of 2-Fi(a, b; c; z) with respect to, e.g. the parameter b, gives (similar expressions are obtained for differentiation with 
respect to a, c): 



2 F^ m (a,b;c; Z ) = £ 



(a) k (b) k ^(b + k) z k 



k=0 



k\ 



-*(6) 2 Fi(a,6; c;z), 



(B.5) 



where ^(fc) is the digamma function and the Pochhammer symbol (a)k = T(a + k)/T(a) (see, for instance, [2J|). With 
the help of formula Eq. (|B.5|) . we get: 

2 F 1 (-£, 2 + e; 1 - e; z) = 1 + e ln(l - z) + 0(e). (B.6) 

Inserting this back into Eq. (|B.4[) . we can write down the result for the integral A m (to first order in powers of x): 



A m {x,y) 



x ^o (x + y + m 2 ) 



2\-l-e 




hi 



to + y J m + y 



+ 0(x 2 )+0(e) 



(B.7) 
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which agrees explicitly with the corresponding expansion of the result given in Eq. (]6 .36[) . 

We now turn to the integral AL, given by Eq. (|6.2p . The parametric form has the expression: 



-4! 



(x + y + m ) 



2\-l-e 



r(l + e) 



(47T) 



2-e 



d0 



1-/3 







(a + /3) 3 / 2 



a/3 



(a + /3) + y + m 2 ) 



(3(1 - (3)y + (3m 2 
x + y + m 2 



(B., 



Calculations similar to the integral A m bring us to the following result (to first order in x): 
x ^o , (x + y + m 2 )- 1 ' - I ^ ' ' 1 



A m { x iV) = h 



(4tt) 2 



2x 

3y 



m 

y 




2 In 



+ 0(x 2 )+0{e) } , (B.9) 



again in agreement with the expansion of the result given in Eq. (j6.48|) . 

Another useful check comes from the study of the mass differential equation. With I n given by Eq. (|6.3|) . the 
derivative with respect to the mass reads: 



Qjn 

am 



uj 2 [(k — uj) 2 + m 2 ]uj 2 



(k — ui) 2 + m 2 



From the relations Eq. (|6.4p . Eq. (|6.5p and Eq. (|B. 10|) we get the following relation: 

Qjn Qjn Qjn 

afc4 afcfc om 



(B.10) 



(B.ll) 



Using the same procedures as in the text, we can then derive a differential equation for the integral in terms of the 
mass: 



, dl n 
dm 2 



(d + ra-4)- 



k 2 



[(k — cij) 2 + m 2 ] 2 a; 2 



(B.12) 



Starting with the case n = where 1° = A m , we see that by inserting the solution, Eq. (|6.36p . we have that in the 
limit e — > 



> dA m 
dm 2 



+ (l + 2e)- 



:A m — 



m 2 (x + y + m 2 ) 



(4tt) 2 



:ln 



'1+221 + 1 



■0(e). 



1 + ^-1 



In terms of Schwinger parameters, the explicit integral of Eq. (|B.12|) reads: 



r(i + e ) f 1 



k 2 + m 2 J [(k — lu) 2 + m 2 ] 2 uj 2 x + y + m 2 (Att) 2 6 J 
and for m 2 ^ indeed 



da (1 — a) 



-l/2-e 



(m 2 + ay) 



-l-E 



(B.13) 



(B.14) 



m (x + y + m ) 



2\-2-e 



k 2 +m 2 J [(k - uj) 2 + ' 



(47T) 



2-e 



:ln 



'1 + 121 + 1 



000 



'1 + ^-1 



(B.15) 



showing that the mass differential equation is satisfied. For m 2 = 0, the right-hand side of Eq. (|B. 13|) vanishes as 
m 2 lnm 2 , whereas the parametric form of the integral in Eq. (|B. 14|) goes like m 2 /e. The integral of Eq. (|B. 14|) does 
contain an ambiguity in the ordering of the limits m 2 — > and e — * 0, but this problem is not of importance because 
of the overall factor m 2 in the differential equation. In fact, since the solution of the mass differential equation is in 
principle formally derived as the integral over m 2 and m 2 = is the only the limit of this integral, the ambiguity 
encountered may be regarded as an integrable singularity and presents no problem. 
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Turning now to the case n = 1 where I 1 = A^, we first extract the overall k 4 factor as before by defining A A = k 4 A m 
such that the differential equation is 

= 2c m2 ~ Xm _ m2 f duj (B.16) 

dm 2 k 2 + to 2 m k 2 + to 2 J [(k — uj) 2 + m 2 ] 2 u) 2 

Notice that in the integral term we have used the identities 

du uj 4 f duj(k 4 -u; 4 ) f duj m-iTi 



[(fc- W ) 2 + m 2 ] 2 J 2 7 [ W 2 +TO 2 ]2 ^_^ 2 J {(k~L0) 2 + m 2 } 2 L0 2 

Now, for m 2 ^ 0, the integral term of Eq. (|B.16j) is finite as e — > 0; however, the m 2 = limit is again ambiguous but 
as above this can be regarded as an integrable singularity. Also, when m 2 — 0, A is known to be e finite (it is the 
massless integral considered in Ref. [la]). This means that as e — > we have the simple integral expression 

™ 2 !% = " 277^2 In 1 ±= I (B.isi 

dm 2 x + y + m 2 (Att) 2 - 6 i i 

Knowing the solution, Eq. (|6.48|) . it suffices to show that when m 2 = the original massless integral from Ref. [15( is 
reproduced and that the derivative of the massive solution satisfies the above. Both of these steps are straightforward. 
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